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Abstract——This paper considers the combined effect of longitudinal diffusion and external 
mass transfer resistance on the performance of fixed bed processes with a favourable non-linear 
isotherm. It is shown that, whereas the linear additive principle is obeyed when the isotherms 
are“ almost linear,” for non-linear systems the two mechanisms are definitely not linearly additive. 
A simple method is proposed for estimating this combined effect with good accuracy. 


Résumé— Cet article traite de Veffet combiné de la diffusion longitudinale et de la résistance 
externe au transfert de masse sur la performance des procédés utilisant un lit fixe avec isotherme 
non linéaire, Ce systéme nobeit plus au principe additif linéaire quand les isothermes ne sont pas 
linéaires ou “ presque linéaires.” L auteur propose une méthode simple pour prévoir leffet 
combiné avec une bonne précision. 


Zusammenfassung —In dieser Arbeit wird die Cberlagerung ciner longitudinalen Diffusion 
und eines dusseren Stoffiibergangswiderstandes auf die Ausbeute cines Festbettprozesses besonders 
mit nichtlinearer Isotherme crortert. Es wird gezeigt, dass bei * nahezu ” linearer Isotherme das 
Prinzip der linearen Additivitat erfillt wird, wahrend sich bei nichtlinearen Systemen dic bieden 
Mechanismen sicherlich nicht linear addieren. Eine einfache Methode wird vorgeschlagen, die 
diesen, Effekt aus der Cherlagerung mit guter Genauigkeit abzuschiitzen gestattet. 


INTRODUCTION which this chromatographic band travels along 
the column, and its relative width. It has 
already been established that the speed is a 
function primarily of the isotherm and can 


Ir nas long been recognized that external resis- 
tance to mass transfer, longitudinal diffusion, 
and the form of the equilibrium isotherm play 
an important role in determining the shape of 
the breakthrough curve for the standard fixed 
bed operations. It is understandable therefore 
that a considerable effort has been spent in an 
attempt to determine to what extent these 
factors affect the performance of fixed bed 
processes and under what conditions they must 
be taken into account in design calculations. 
The more noteworthy contributions on this 
subject have been recently reviewed in the 
comprehensive article by VermMeuLen [1]. 

Two general characteristics of breakthrough 
curves are of particular interest; the speed with 
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therefore be obtained with good accuracy from 
the equilibrium theory. On the other hand, for 
isotherms of the so-called favourable type (see 
Fig. 1) the relative width is not seriously affected 
by the form of the equilibrium curve, but is 
primarily a function of the finite rate of mass 
transfer and of longitudinal diffusion, if other 
usually secondary effects, intraparticle diffusion 
for example, are neglected. 

Earlier studies, far too numerous to be men- 
tioned here in detail, were concerned with the 
problem of determining to what an extent these 
two mechanisms, taken separately, were res- 
ponsible for the band spread often observed 
in chromatographic and ion-exchange columns. 
It appears, however, that for some processes the 
effects of a finite rate of mass transfer and longi- 
tudinal mixing are comparable, which therefore 
raises the interesting question of how these two 
mechanisms combine. Are the effects linearly 
additive, as has been sometimes postulated, or is 
the physical phenomenon far too complicated for 
such a simple solution ” 

This problem was recently investigated by 
et al. [2,3] who restricted them- 
selves to a linear isotherm. By simplifying a 
rather general mathematical solution first obtained 
by Amunpson and Laprpus [4] these authors were 
able to prove that indeed the two mechanics of 
longitudinal diffusion and gas-phase mass transfer 
resistance are exactly additive and that quantities 
such as the H.E.T.P. have the linear additive 
property. The same result was also obtained from 
a statistical model [2] under the assumption that 
the two mechanisms are statistically independent. 

These conclusions reached by KiLuNKENBERG 
et al. apply of course only to dilute solutions for 
which the isotherms may be assumed linear, 
and one would certainly not expect that for 
non-linear systems the two mechanisms would 
act independently. The purpose of this note 
therefore is to pursue this study somewhat 
further and consider the combined effects of 
longitudinal diffusion and external mass transfer 
resistance on the performance of fixed bed 
processes for non-linear isotherms of the favour- 
able type. As will be shown, this problem admits 
of a relatively simple mathematical solution. 


ror Favourasie Basic EQuatrions 
AND THEIR LsoTHERMS 
According to the well-known basic equations 
for a one-dimensional fixed bed process, 
a 


10-09% 
Pom! 


vt 
if isothermal conditions are assumed, and if 
effects other than those caused by longitudinal 
mixing and by a finite rate of external mass 
transfer are neglected. The symbols have the 
following meaning : 


¢ = concentration in moles per unit volume 

of absorbate in the moving stream. 
- concentration in moles per unit volume 

of the stationary phase. 
the combined molecular and eddy long- 
itudinal diffusion coefficients. 
the fractional void volume of the bed. 
the superficial velocity of the fluid 
through the bed. 
a mass transfer coeflicient with reci- 
procal time units. 
the equilibrium value of ¢ corresponding 
to a given stationary phase concen- 
tration q. 

X = distance along the axis of the bed. 

t — the time 


and the boundary conditions will be taken, for 
simplicity, as 


= ¢gandg — Qfor X — where c* (Q) = eg. 
+ + ow. 
¢,H(— X)jatt—0. 


where H(z) is the Heaviside step function such 
that H(z) = 0 for: < O and H(z) = 1 forz > 0. 

As might be expected, a complete solution of 
equations (1) and (2) has not yet been discovered, 
except for the special case of a linear isotherm 
4). It was first noticed by and 
Exepaut [5] however, that in the absence of 
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longitudinal diffusion (D = 0) and for an iso- 
therm of the so-called favourable type (Fig. 1) 
where 


0 for allq > (3) 
the breakthrough curve approaches a constant 
pattern as it moves down the column. A similar 
result was recently reached by Licurroor [6] 
for the case of an infinite rate of mass transfer 
but a non-zero longitudinal diffusion coeflicient. 
Of course, the constant pattern is only an asym- 
ptotic profile which is strictly speaking attained 
only at distances far removed from the entrance 
of the column. Yet it has been established by 
Laripus and Rosen [7] that, as is often the case 
for analogous physical problems, such a constant 
pattern is quickly reached inside the bed and may 
therefore be used in design calculations with good 
accuracy, unless the column is very short. 

It is easy to show now that ¢ and g become 
asymptotically functions of a single variable, if 
the isotherm is of the favourable type. Thus, if 


— Ucgt 

~ eg + (1 Q] Dlecg + (1 — 
then, asymptotically, equations (1) and (2) reduce 
respectively to 


(4) 


du dv 
dz dz 


which can immediately be integrated into 


u*(v)—u 


u*(v)—u 
u—v 


ve! and ut (v) = 
Q 0 


These equations also must satisfy the boundary 


(8) 


w=v=0 at 


and by definition, 
u=}jatz=0 (9) 


We see then that the dimensionless group <. 
defined as 


—«€)Q Dk 


plays indeed a very significant role in the per- 
formance of fixed bed operations and measures 
in a sense the relative importance of longitudinal 
diffusion and external mass transfer resistance in 
causing band spreading. 

Again, a closed form analytic solution of 
equations (5) and (7) cannot, unfortunately, 
be obtained for an arbitrary isotherm u* (v) 
unless x +0 or 

For x u u* (v), and therefore 


(10) 


where u* (v*) =u (11) 


as obtained by Licurroor [6]. On the other 
hand, if « u +v, and 


- du 
u (12) 


However, for intermediate values of x a 
numerical integration of the equations is indicated. 
Nevertheless, it is possible to arrive at an approxi- 
mate solution, perfectly adequate for most 
applications, by elementary methods. 


Approximate solutions 


Let 
(v) = rv + g(v) (13) 


where r is equal to the slope of the isotherm at 
v= 0, 
and v Sau + v, (u) (14) 


so that g (0) = g’ (0) = v, (0) = v,' (0) = 0. Upon 
substitution in equations (5) and (7) one readily 
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Therefore 
dv 
with 
finds that 
conditions 
e*=v=1 at 
and that 
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=(a— I)u + v, (4) (16) 
Therefore, since by definition u = 4 when z = 0, 
o5 
(uy) (a—1) Qu 


a— (a 1) u? ur, du (17) 


In particular, if u +0 


fing. 4| (18) 
x} | 


A=- | du (19) 
0 

Of course, an entirely similar solution can be 
obtained for u +1. However, since one is in 
general interested primarily in the distance z 
between the midpoint of the breakthrough curve 
(wu = 1 2) and a small value of u, u = 0.05 say, 
equation (15) is usually quite adequate. The 
constant A can be determined accurately only 
by either a numerical solution of equations (7) 
and (19), or by a series expansion. Fortunately, 
though, a fairly reliable estimate of A can be 
arrived at by considering its asymptotic forms 
for some typical isotherms and then interpolating 
between the two extremes. 


Tue Paramerer A vor Constanr Excuance 
Facror 
For a constant exchange factor, the equilibrium 
isotherm becomes 
wv 
u* (v) = Ferre (20) 
where +, <1 if the isotherm is_ favourable. 
Equations (5) and (7) can then immediately be 
solved for the extremes « + 0 and x - «. Thus, 


(a) For «+0 


(2u)’ 


l-r 


and therefore 


as x +0 (21) 


(b) For z+ @ 


in [2(1 — 


so that 
A > (22) 


For intermediate values of x, a simple inter- 
polation formula can readily be obtained. One 
notes from equation (15) that the two mechanisms 
have a comparable effect on the overall solution 
where r ~ , and it appears reasonable therefore 
to expect that 

In 2 1 

which reduces to equations (21) and (22) as 
x and « respectively. It should be 
remarked here that formulae similar to that 
given by equation (23) are not infrequently used 
in physical problems to interpolate between two 
exact asymptotic expressions, and that the inter- 
polation formula chosen here is therefore not as 
arbitrary as it might appear at first glance. 
It is seen from Table 1 moreover that the agree- 
ment between an exact numerical solution for 
r— 1 2 and equation (23) is indeed excellent. 


Table 1. A comparison of the exact values of A for 
r — § and those calculated from equation (23). 


A (exact) A from eqn. (23) 
0 139 1:39 
0-87 O-87 
4 | O54 0-56 
1 Ow O46 
2 0-38 O-41 


As a point of interest we might mention here 
that as r +1, an exact solution may be arrived 
at without difficulty. Thus, in the limit as r —- 1. 


(l—r)z_ 2(1—u) 


+o 
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for all x, which clearly shows that the two 
mechanisms of longitudinal diffusion and external 
mass transfer resistance are indeed linearly 
additive for “almost linear systems” only. A 
similar conclusions would of course be reached 
for a system with any “* almost linear ” isotherm. 
Naturally, a similar approach can also be used 
successfully for other isotherms. In general the 
two extreme cases (x +0, or x +» ©) can be 
solved directly without too much effort, and the 
interpolation formula, equation (23), would be 
expected to hold with a fair degree of accuracy. 
It is believed for example that such an approach 
could be used to investigate the effect of longi- 
tudinal diffusion on the behaviour of a chromato- 
graphic column with the rather commonly 
employed Thomas-type kinetic law [8, 9}. 


Tue Breakrurovcnu (r = 0) 


Equation (23) appears to be accurate for 
1/4<r<1. As however the irreversible case, 
r +0, is approached, a different solution is 
required, Fortunately though the equations can 
be solved exactly for this special but rather impor- 
tant problem. Thus, if 


u*(v) = 0 for O <1 
and 


u* = u for v 1 


it can be shown by elementary means that : 


a) 


2A 
z=In2(1~— u) for In (26) 
and 
» 
2A 
for z >In wer} (26a) 


(b) If on the other hand A < 1, then 


On the combined effect of longitudinal diffusion and external mass transfer resistance in fixed bed operations 


A+1{ A+ 
wt 
for = <In (27) 
and 
1 1 2A 
for z > In (27a) 


We easily see therefore from equations (18), 
(26a) and (27a) that, when r = 0, 


2 
In 


A= + 


for A>1 (28) 


A =0 for A <1 


It might be worth remarking here that an analytic 
solution similar to that shown above, can also 
be readily obtained for an isotherm of the form 


u* (v) = rv for v <1 
= for v = 1. 


This solution will however not be reported here 
since such isotherms do not appear to be very 
realistic. 


CONCLUSIONS 

A simple procedure has been presented for 
estimating with good accuracy the combined 
effect of longitudinal diffusion and external mass 
transfer resistance on the performance of fixed 
bed processes with a non-linear but favourable 
isotherm. This analysis shows that the linear 
additive principle is obeyed when the isotherm 
is “almost linear” but that for non-linear 
systems the two mechanisms are definitely not 
linearly additive and must be taken into account 
by means of equations (18) and (23) for 
}<r<1, and equation (27) for r +0. It 
appears, therefore, that conclusions arrived at by 
a study of linear problems about the way in 
which the different mechanisms combine in fixed 
bed operations should not be automatically 
extended to non-linear processes, without at 
least some additional work. 
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Study on cooler condensers for gas—-vapour mixtures 
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Kyoto University, Kyoto, Japan 


(First received 16 November 1959 ; in revised form 1 April 1960) 


Abstract—The graphical method for the design of cooler condensers by Mizusutna et al. presented 


in a previous paper is revised. 


Furthermore gas- vapour mixtures, such as water air, benzene -air, toluene—air and methanol- 
air were cooled in a vertical single tube condenser, and the experimental data were found to be in 
good agreement with the results of calculation by the method of this paper. 


Résumé—Les auteurs révisent la méthode graphique de calcul de condenseurs de Mizusuina 


et Col, déja publiée. 


De plus des mélanges gaz-vapeur tels que eau—air, benzene air, toluéne—air et méthanol-—air 
sont refroidis dans un condenseur a tube vertical unique, et les données expérimentales sont en 
bon accord avec les résutlats des calculs obtenus par la méthode donnée. 


Zusammenfassung 


Die graphische Methode zur Planung von Kiihler-Kondensatoren, die von 


MizusHina u.a. in einer friiheren Arbeit veréffentlicht worden war, wurde iiberarbeitet. 
Weitere Gas-Dampf-Gemische, wie Wasser Luft, Benzol-Luft, Toluol-Luft und Methanol— 

Luft wurden in einem vertikalen Einrohrkondensator gekiihlt. Die experimentellen Ergebnisse 

stimmen mit den nach der Methode dieser Arbeit errechneten Ergebnissen gut iiberein. 


INTRODUCTION 


Since the stepwise method for designing cooler 
condensers was published by Coisern and 
HovuGen [1] various design methods have been 
proposed by Bras [2, 3, 4], Crise and Nevson [5] 
and Hutpen [6]. These papers described methods 
of calculation without experimental data. Cairns 
[7] presented the experimental results of the 
water-air and chloroform air systems, and Surru 
and Ronson [8] the experimental data of the 
benzene-air system. In both papers the experi- 
mental results were compared with those of 
and Hovucen’s method. 

Recently, Mizuswuina et al. [9] published a 
graphical design method for cooler condensers, 
which was applicable to any gas- vapour mixture 
system. In this previous paper even in the case 
of mixtures having high vapour concentration, 
it was assumed that h, ky — C, for water air 
system and that pyy (P — p;) — 1 for mixtures 
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other than the water-air system. These assump- 
tions however sometimes cause considerable error 
in estimating heat transfer surface area. 

In this paper a more accurate method involving 
some additional procedures will be developed, and 
the results of the calculation will be shown to be 
in good agreement with the experimental data 
obtained for four systems, i.e. water-air, benzene— 
air, toluene—air and methanol-air in a single tube 
cooler condenser. 


EXPERIMENTAL APPARATUS AND Procepure 


The apparatus consists principally of a boiler 
for the liquids, a vertical single tube condenser, a 
condensate receiver, an air compressor and a 
cooling water pump. 

Air was pumped by a } h.p. compressor through 
a CaCl, column for drying and heated electrically 
before being mixed with the vapour. Its flow rate 
was controlled by a valve and measured with an 
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orifice meter. Various kinds of liquids were 
evaporated in the boiler and mixed into the air 
flow. The boiler had an internal capacity of 30 em 
diameter »« 40 em height, being heated with two 
2 kW immersion heaters, The electric power input 
was controlled by a variac, and measured by an 
ammeter and a voltmeter. The vapour rate was 
measured by the power input. Power input vs. 
the rate of evaporation of liquid was calibrated 
beforehand. The heat loss from the boiler was 
obtained as 1-59 Mf keal hr, where Af C was the 
temperature difference between boiling liquid 
temperature and the surrounding air tempera- 
ture. In the case of the high vapour content 
water air system an auxiliary boiler was used to 
produce the additional vapour. The gas vapour 
mixture was led to the condenser through a 
} in. steel pipe. 

The condenser was a vertical double-tube type. 
The gas vapour mixture flowed downwards inside 
the inner tube, and the cooling water tlowed 
countercurrently in the annulus. The inner tube 
was made of copper and of 29 mm o.d., 23 mm 
id. and the length of the tube exposed to the 
cooling water was 180 em, ic. the effective internal 
surface 
calming section was provided before the condenser 


area was long 
and a 20 em long section was extended below the 
hottem of the The external tube 
forming the water jacket was Esron 
(polyvinylchloride resin) pipe (i.d. 50 mm). 
The condensate was received in the condensate 


condenser. 
2 in. 


receiver, which was attached to the bottom of the 
copper tube. The rate of 
determined directly from the weight of the con- 
densate collected during a known time. 

The cooling water was pumped by a | hip. 
turbine pump and metered with an oritice meter. 

Calibrated copper constantan 
were located at the inlet and outlet of the cooling 
water, the air line before entering the vapour 
line, the inlet and outlet of the gas vapour 
mixture and also at ten points located at equal 
intervals on the axis of the inner tube of the 
condenser. Thermocouples were also set on the 
wall of the inner tube and midway between the 
annulus walls at heights corresponding to those 
on the axis. Thus the change in the temperature 


condensation was 


thermocouples 


of the gas-vapour mixture, the pipe wall and the 
cooling water were measured. The electrical 
potential was measured by a potentiometer 
within an aceuracy of 1 xV which corresponds to 
about 0-03 C. The setting of the thermocouples 
is shown in Fig. 1. 


Stainless stee! 
protecting tube 
Copper protecting 9 


Stainless stee! 
protecting tube 


Fic. 1. Details of setting of thermocouples. 


boiler was insulated with Sem thick 


The 


diatomaceous earth and all air lines and gas 
vapour mixture lines were lagged with 2 em thick 


glass wool, 

In making each run the boiler was filled with 
liquid, the power input to the boiler was set at 
the required vaporization rate and air was intro- 
duced into the vapour line. 

Sullicient stabilizing time was allowed for each 
run, At steady state conditions, the temperature 
at each point, the power input to the boiler and 
the rates of air, vapour, condensate and cooling 
water were measured, The results of the experi- 
ments are listed in Table 1. 


Mernop or 

Since the publication of the previous paper [9] 
on the design of cooler condensers for gas) vapour 
mixtures the authors have calculated the heat 
transfer surface area for many different cases, 
and it was found that the assumptions that 
h, k,’ ~ ©, for mixtures of water vapour and air 
and — 1 for mixtures other than 
the water air system sometimes caused consider- 
able error in estimating the surface area. In 
order to improve the accuracy of the graphical 
method the following additional calculation 
procedures were developed. 


_fube Esron pipe 
ome Me ope’ 
pe 
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Basic equations 

For the sake of convenience the important basic 
equations presented in the previous paper [9] 
will be noted again here. The equations are 
numbered to agree the numbering in the previous 
paper [9), 

The ratio of the heat and mass transfer co- 
ellicient ts 


(3) 


(4) 
G (Pr). 
The enthalpy balance is 
dA (6) 


Gy di, = 


aC, 


where 7’ is the * modified enthalpy ~ and is detined 


as 


Therefore 


Deriving 


where 


equation becomes 


hn 


A heat balance yields the slope of the operating 
line 
(9) 
The slope of the path line of the gas- vapour 
mixture is 


and the slope of the tie line is 


_ _ Cy 
hn 


The difference between the enthalpy and the 
modified enthalpy is given by 


(x —1)C,4. 


i,’ — i, 


 — 


ig (23) 
The additional procedures 

(1) Assuming that the mixture is saturated at 
the exit, caleulate the values of 8 at the inlet 
and exit, 

(2) Taking as a first approximation x — 1, ie. 
h, ky Cy. construct tie lines of slope 
— Bhy Cy hy through the points A and Z on the 
operating line and locate A” (t,o, and 
(tig on the saturation curve as shown 
in Fig. 2. 


Alig, .ty Hy) 


ty) 


ta 


Fic. 2. Enthalpy and modified enthalpy— 
temperature chart. 


(3) Calculate the value of a» and ag» by 
equation (4) using values of p, and py which 
correspond to t; and ty respectively. 

(4) Take the arithmetic mean value of x, 
and az for «. Usually, the value of x is less 
than one. 

(5) Using the value of « obtained, follow the 


procedure in the previous paper [9]. 


h, 
k, 
where 
) Seyi 
VOL. 
13 
96 
aC, t, + Ag Hl, lai 
2G | a 
(19) 
C BC, : 
B=(1+ H,)" (20) 
= 
di, Vi,’ 
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Tur Comparison or Resuts or 
CALCULATION AND EXPFRIMENT 


The procedure mentioned above was applied to 
all of the experimental runs. An example of such 
a calculation is given in the following section. 


Run No. 1 (water-air system) 
At the gas inlet: 


=~ 510 kg kg of dry air 
~ O-418 keal kg 
O95 
47200 kg hr 
20970 


Since the preliminary experiments on the heat 
transfer coefficient of this apparatus gave 


j = 0-0822 Re ®?, 
j = 00044 


Substituting these values into equation (12) 
the heat transfer coetlicient of the gas film h, is 


calculated as 


h, = (O-0044) (0-418) (47200) (0-95)t — 
89 keal m? hr C 


The value of C, obtained from Fig. 1 in the 
previous paper [9] is 


C, = 2-55 keal kg of dry air C 


Since the gas flows inside the tube n — 0-2 in 
Fig. 1 in the reference [9]. Thus the value of 8 
is obtained as, 

1-436 


From equation (19) 


hyn Bh, (1-486) (89) 

Since the heat transfer coeflicient of the cooling 
water was found to be 6900 keal m?® hr C in the 
preliminary experiments, and since the vapour 
film coeflicient is calculated as 3260 keal_ m? hr °C 
from Nusselt’s equation, the value of the combined 
conductance is 


hy = 2210 keal/m? hr °C 


Since the gas condition at the inlet is 
= 97-9 C, Hy = 5-10 kg kg of dry air, the 


enthalpy of the gas-vapour mixture at this point 
is calculated as 


i, — Cyt, + Ag H, = (2-55) (97-9) + (595) (5-10) 
~ $285 keal kg of dry air 


Assuming the mixture at the outlet is saturated 
at 698 


H, = 027 kg kg of dry air 

Cy 0361 keal kg of dry air 
i, = 187 keal kg of dry air 

The temperatures of the cooling water are (see 


Table 1), 


at the gas inlet ¢, = 14-8 °C 
at the gas outlet ¢, 85 C 


Thus the operating line can be plotted on the 
enthalpy temperature chart as shown in Fig. 2. 

Assuming h, k,’ — C,, the slope of the tie line 
is calculated from equation (22'). 

At the gas inlet 


Cp hy — — (1-436) (2210) (0-0200) 
63-5, 
At the gas outlet 
Bho Cp hy = — (1-049) (2210) (0-0200) — 
~ 
From points A and Z on the operating line tie 
lines with these slopes are constructed and the 
intersections A’ and Z'” give the interface 
temperatures, 
tg = 64-2 °C 
tiger = 124°C 
Hence, 181 mm Hg, 
= 108 mm Hg. 
Since at the gas inlet p,, — 677 mm Hg, 


and at the gas outlet p,z — 232 mm Hg, 


then 


(Pum)a 
(760 — 181) — (760 — 677) 


— = 255mm Hg 
2-308 log (760 — 181), (760 — 677) 


— Pi/a 579 
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Similarly 
(Pam)z ~ 631 mm Hy 


(57 Pam 0-843. 


The Schmidt numbers are 0-61 at the gas inlet 
and 0-62 at the gas outlet, and the Prandtl 
numbers are 0-95 at the gas inlet and 0-79 at the 
gas outlet. 
From equation (5) 


(0 | 


~ (0-843) (— = 4) 


(0-35 + 0-75) 


2 


Hence x 


The relation between (i, — 7,') and ¢ at various 


values of H, is computed by equation (23) using 
x — 055 and plotted in Fig. 2. The saturation 


curve for 7’ is also constructed, 


i004 


Woter-air 


Fic. 3. Comparison of the results of calculation and 
experiment for the gas temperature change. 


Water Air 
Run No. 
~ Run No. 
© Run No. 
@ Run No. 


Calculation by the method of 
this paper 


Experimental data. 


The slope of the tie line is given by equation 
(22’) as 


zBhyCy hy (0.55) (2210) (0-0200) 8 
24-3 8 

At point 4 of Fig. 2 
il 


1-436, i, 


5-10, 


3170 is obtained from Fig. 2. 
xBhyCy hy 


and iy 


Hence (243) (1-486) - 


— BES 


straight line from point a 3170, 
t; ~ 148 C) with a slope of — 348 gives the 
interface conditions at the gas inlet, Le. i;° — 645 

The slope of the path line is obtained from 
equation (10°) 

di 1 3170 — 645 


- = 430 
dt, 0-55 97-9 — 87-2 


The intersection of the straight line from the 


© 002 004 006 008 00 O14 
A, m 


Fic. 4. Comparison of the results of calculation and 
experiment for the gas temperature change. 


Benzene— Air 
- Run No. 
- Run No. 
© Run No. 1: 
@ Run No. 


Calculation by the method of 
this paper. 


Experimental data 
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point A’ (i, = 3285, t, = 97-9) with a slope of 
430 and the line of constant temperature 92 °C 
gives a point B’ (i, — 1183, t, — 92-0). 

These steps are repeated for the remaining 
points. The results of point to point calculation 
are shown in Table 3(a). The value of the integra- 
tion in equation (21’) is 3-79. The surface area is 
therefore 


(0°55) (0-0200) (3-22) (3-79) — 0-134 m2 


Since the surface area can also be calculated by 


— 
LC, at 
ho | —t 
— (1445) (19) 194) — 0-197 mt 
2210 


equation (21°) may be simpler to use than 
equation (21'). 

Comparison of the calculated results with the 
value of 0-130 m? used in the experiment indicates 
that the agreement is good, 


| 
0 002 004 006 008 010 O14 


A, me 


Fig. 5. Comparison of the results of calculation and 
experiment for the gas temperature change. 


Toluene— Air 
—— Run No, 21 
Run No. 23 
© Run No. 21 
@ Run No, 23 


Calculation by the method of 
this paper. 


Experimental data. 


The same method of calculation was applied 
to every run of the four systems, i.e. water—air, 
benzene air, toluene-air and methanol-air, the 
results being tabulated in Table 2. The results of 
the calculation coincide well with the experimental 
data. The change in the measured gas temperature 
in the condenser is compared with the results of 
calculation in Figs. 3-6. Again there is good 
agreement between the results of calculation and 
the experiments. 


Discussion 


The comparison of the calculation results of 
the method described above for the water air 
system (Run No. 1) with those of Coupern and 
Hovcen [1] and Hunpen [6] is shown in Table 
3(a). The surface areas obtained by these three 
methods were 0-134 m*, 0-136 and 0-128 
respectively, and these results showed good 
agreement with each other. The results of cal- 
culation for the benzene air system (Run No. 11) 
are also tabulated in Table 3(b). The agreement 


is also good. 


— 

Methanol! - air | 
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Fic. 6. Comparison of the results of calculation and 
experiment for the gas temperature change. 


Methanol Air 
—— Run No. 31 
Run No. 33 
oO Run No. 31 
@ Run No. 33 


Chem, Engng. Sci. Vol. 13, No. 1. August 1960. 


Calculation by the method of 
this paper. 


Experimental data. 
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In the calculation of hy the vapour film co- 
ellicient is obtained from Nusselt’s equation. For 
the water-air system where the vapour film 
coetlicient is very high an average value of the 
coetlicient throughout the condenser may be 
used. However, for systems other than water—air, 
where the value of the vapour film coeflicient is 
relatively small, the point value should be taken. 
For the methanol-air system (Run No. 31) the 
value of hy varied from 6600 at the gas inlet to 
800 at the gas outlet, and the surface area obtained 
from the calculation using the point value of the 
vapour film coeflicient was 0-121 m*. The surface 
area obtained from the calculation with a constant 
average value hy — 1060 was 0-122 m*, the two 
results for the surface area being nearly equal. 
Therefore, as shown in Table 2, the values for h, 
were taken as constant throughout the condenser 
for all systems. 

If the amount of heat to be transferred is 
calculated from the enthalpy change of the gas, 
it includes the heat necessary to cool the conden- 
sate to 0'-C, Actually, the temperature of the 
condensate lies between t, and t,. From this point 
of view, the graphical method using the enthalpy 
chart of the gas gives a conservative estimate of 
the heat transfer surface area. 


Cs humid heat keal/kg of dry gas °C 
F factor defined in reference [6] 

Gy flow rate of stagnant gas kg of dry gas/hr 
G mass velocity of gas vapour mixture — kg/m? hr 
H humidity kg/kg of dry gas 

h coefficient of heat transfer keal /m? hr °C 
hy fictitious coefficient of heat transfer when 
only stagnant gas flows keal/m? hr °C 
kg combined conductance other than the 
gas film keal m?® hr °C 
i enthalpy of gas vapour mixture keal kg of dry gas 
i = modified enthalpy of gas vapour mixture 
keal kg of dry gas 
jj factor in heat transfer 
k, transfer coefficient 
kg ‘m* hr (kg kg of dry gas) 
ky mass transfer coeflicient kg m* hr mm Hg 


tlow rate of cooling water kg hr 
P total pressure in condenser mm Hg 
p partial vapour pressure mm Hg 
Pym average partial pressure of stagnant gas 
in gas film mm Hg 
temperature c 
constant defined by equation (4) 
variable defined by equation (20) 
latent heat of vapourization at 0° C keal kg 


hypothetical temperature defined in reference [6] 
Prandtl number 
Reynolds number 
Schmidt number 
Subscripts 
B stagnant gas 


Fed wa 


Noration 2 gas-vapour mixture or gas film 
A area of transfer surface i interface 
C specific heat keal kg °C {cooling water 
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Optimum temperature gradients in tubular reactors 


R. Aris 
Department of Chemical Engineering, University of Minnesota, Minneapolis 14, Minnesota. 


(Received 1 December 1959) 


Abstract—The method of dynamic programming is used to solve the general problem of finding 
the temperature gradient in a tubular reactor, which will maximize some function expressing 
the profit made by the reaction. As an introductory example the reaction A > B->C when 
the yield of B is to be maximized is considered. This is followed by a detailed analysis of this 
reaction for first order kinetics, which confirms and extends the work of AmuNpson and BiLous. 
The solution of the general problem for a number of simultaneous reactions is then presented 
and some extensions discussed. 


Résumé——La méthode de programmation dynamique est utilisée pour résoudre le probléme 
général de la détermination du gradient de température dans un réacteur tubulaire, qui portera 
au maximum la fonction exprimant le progrés d'une réaction. L’exemple préliminaire est celui 
de la réaction A ->» B -» C quand le rendement de B est poussé au maximum. Suit une analyse 
détaillée de cette réaction pour une cinétique d’ordre 1 qui confirme et étend le travail d’- AMUNDSON 
Bitovs. La solution du probléme général pour un certain nombre de réactions simultanées est 
ensuite présenté et son extension discutée. 


Zusammenfassung -Die Methode der dynamischen Programmierung wird benutzt, um das 
allgemeine Problem zu Wisen, den Temperaturgradienten in cinem Rohrreaktor zu finden, bei 
dem eine Funktion den Reaktionsfortschritt ausdriickt, die den Maximalwert annimmt. Als 
einfiihrendes Beispicl wird die Reaktion AB —> C betrachtet, bei der die Ausbeute von B ein 
Maximum erreichen soll. Es folgt cine ins EKinzelne gehende Analyse dieser Reaktion unter der 
Voraussetzung, dass die Teilschritte nach der ersten Ordnung verlaufen, Damit wird die Arbeit 
von AmuNDsON und BrLous bestitigt und erweitert. Im Anschluss daran wird die Lésung des 
allgemeinen Problems fiir eine Reihe von Simultanreaktionen angegeben, und einige Erweiterungen 
werden diskutiert. 


1. InrropuctTION 


Ir nas long been known that for a single reaction 
taking place in a tubular reactor the best yield is 
obtained if the temperature is so adjusted as to 
make the reaction rate as large as possible at 
each point of the reactor [1, 2, 3, 4]. In the case 
of two or more competing reactions this no 
longer true, as AmMuNnpson has shown [5]. The 
problem of two consecutive reactions A + B + C 
has been extensively discussed in the work of 
Bitous and Amunpson [4] and they derive an 
equation to be solved simultaneously with reaction 
equations to give the optimum temperature 
profile. They use two methods, the first based on 
an explicit solution of the reaction equations, 


and the second an implicit method which may be 
applied to the reaction equations themselves. 
Both involve the rather difficult notions of 
functional differentiation [6] and, perhaps because 
of this, their work has not received the attention 
it deserves. 


In the present paper we wish to show that these 
problems may be solved by the methods of 
dynamic programming. These have been 
developed in quite other fields of study but the 
new formalism of the calculus of variations given 
by Be.iMan is extraordinarily apt for just such 
problems as we are concerned with. BeLuMan 
({7] Chap. IX) discusses the maximization of an 
integral 


T 


J (2) = 2) dt (1) 
0 


where the function z(t), 0 <t< T, has to be 

chosen to give the maximum value of J and the 

a, are determined by the differential equations 
dx, 


dt + (0) = (2) 


i= 12m 
and subject to quite general restrictions 
Ry, 2) SO. (3) 


Now if we think of the 2; as concentrations or 
partial pressures, z as the temperature, and f 
as the current holding time along the reactor 
whose total holding time is 7, his formulation is 
just what we need. For if, for example, we wish 
to maximize the yield of the first component, 
putting F = G, gives J = 2,(T) — 2, (0), the 
increase in its concentration. The equations (2) 
are the equations for the rates of reaction and the 
conditions (3) include such restrictions as insist- 
ing that the concentrations be positive and 
bounded (0 < a, < 2,;*) or that the temperature 
be within certain limits (z, <2 < 2*). We shall 
apply BeLuMan’s methods to the general problem 
and give a detailed treatment of the case 
A-+>B-C which Amunpson and BiLous 
studied, incidentally resolving two questions 
that they left open. 


2. Tue Dynamic ProGRAMME FoR A 
Continuous Decision 


The determination of the optimum temperature 
gradient may be regarded as making a continuous 
decision, At each point the temperature has to 
be chosen so as to maximize the eventual yield 
of the process. This choice will depend on the 
concentrations at the point and on the holding 
time of the remainder of the process, and is 
called the optimum policy relative to this state 
and holding time. Now whatever policy may have 
been adopted up to a given point the overall 
optimum cannot possibly be attained unless the 
part of the reactor beyond this point uses the 
optimum policy with respect to the state at that 
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point and the remaining holding time. In fact 
if we vary the policy in the first section and 
always use the optimum policy in the second 
section the maximum over all such variations 
will give the optimum for the whole. But in 
this case we have only to vary the temperature 
in the first section, which can be as short as we 
please. This suggests a way of building up the 
solution and we will illustrate it by a definite 
example. Suflice it to say that the basis of the 
method has been enshrined by Bettman in the 
so-called principle of optimality, that “an 
optimal policy has the property that whatever 
the initial state and initial decision are, the 
remaining decisions must constitute an optimal 
policy with regard to the state resulting from the 
first decision.” ([7], p. 83). This simple statement 
is worth weighing ; its proof is immediate but it is 
pregnant with the most diverse applications. 

Consider the consecutive reactions A + B + C 
with C a waste product and B the product of 
value; as Bitous and AmuNbDson remark ([2] 
p. 88) several reactions may be included under 
this form. If # and y are the concentrations of 
A and B the kinetics of the reactions are given 
by 

dx 


Fe) (4) 


—na(T) F(x) —b(T)G(y) (3) 


where ¢ is the holding time of the reactor up to a 
given point and T7'(t) the temperature there. 
a(T) and 6(T) are the rate constants of the 
reactions. The yield of B over the total holding 
time ¢ = S of the reactor has to be maximized, 
and this is 


~y(o) = | 
0 


s 


= F(x) (6) 


When the optimum temperature policy T (t), 
0 <t < S, has been found the resulting maximum 
will be a function of x. yo, and S. Let 
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Max | F (x) — b(T) G(y)} dt (7) 


then dividing the total holding time interval 
(O, S) into a very short first section (O, s) and 
the remainder (s, 8) and applying the principle of 
optimality 


Yo 8) = 
Max | — S(a(s). y(s).8 (6) 


where in this maximization we have only to 
choose T(t) over the interval (O, s). Now let s 
be very small then 


S Yo. S) = Max oe: — s f(x Yo 8) 


where a suffix 0 denotes values at the inlet 

= Oand the maximization is now over the choice 
of T (s) at s = 0, ie. Ty. In this equation second 
and fifth terms on the right side are unaffected 
by the maximization and the second is equal to 
the left side. Thus using (4) and (5) the equation 
can be rearranged and letting s +0 we have 


This maximum must be sought subject to the 
physical restriction T, < T, < T*, where T, 
and 7* are the minimum and maximum per- 
mitted temperatures. It is a valuable feature 
of dynamic programming that it is not at all 
embrassed by such restrictions; conventional 
calculus of variations is more bashful. 

If the value of 7, is free of restrictions the 
optimum value may be found by differentiating 
equation (9). Thus 7, should satisfy 


0 = {na'(T,) F (x) — (Ty) G (yo) } 
yf of 
(1 + F (x9) = = (10) 


and a partial differential equation of the first 
order for f would result from eliminating Ty, 
between equation (10) and 


¥ = {na (Ty) F (x) — 6(T,) G (yo)} 
F (0) (11) 


Since rate constants such as a or b are generally 
monotonic functions of temperature it is possible 
and advantageous to take one of them as the 
independent variable. Let 6 =<b(a) and 
db da = b'(T) a’ (T) A(a), thus a has to be 
eliminated between 


fs = (naF — b(a)G)(1 + f,) — aF f, (12) 
and 0 = (nF — A(a)G) (1 +f,)—Ff, (18) 
(Here the notation fy = f, = 


Sy = Ve is used for simplicity). 


Suppose that the equation obtained by eli- 
minating a between (12) and (13) is 


H (fy. Sy Yo) = (14) 
for clearly f itself and S do not occur, The 
characteristic equations of this partial differential 
equation are 

diy dy dS dH 
dr dr dr fs 
df, 

dr dr yy’ 
where 7 is a parameter along the characteristic. 
It may be well to notice the importance of the 
last equation of (16) at this point. f, is the partial 
derivative of the optimum profit with respect to 
reactor length S, and the equation shows that 
this does not change along a_ characteristic. 
Naturally we are only interested in the region in 
which f, > 0 for if fg <0 it is unprofitable to 
continue the reaction. 

If the elimination of a between (12) and (13) 
is explicity possible we arrive at an expression 
for H and thence to the characteristic equations 
(15) and (16). However, even if this is not pos- 
sible the characteristic equations can be found 
from the two equations by differentiating them 
partially with respect to f,, f, ete. and eliminating 
the partial derivatives of a. This gives 


(15) 


— 
— 
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-—aF (17) 


dye = nak —bG 


de (18) 


(19) 
=aF'{f,—n(l+f,)} (20) 


(21) 


= bG' (1 + f,) 


These equations are not complete in themselves 
unless (13) can be solved for a and a substituted 
in them, but in this case the elimination would 
be explicitly possible. However, we can find an 
equation for a along the characteristic by differen- 
tiating (13) with respect to + and substituting from 
(17-21). This gives 

P @ (Yo) 


— = — aB(a) 


; (22) 
dr G (Wo 
where 


B(a) = {b(a) —aA(a)} A’ (a). (23) 


It is at first sight rather remarkable that f, and 
J, do not appear in equation (22), for it now allows 
the four equations (17)-(19), (22) to be solved 
simultaneously. We shall show later that this 
always comes about with two simultaneous 
reactions but is not to be expected with more 
than two. 

It is interesting to compare these four equations 
with those of Amunpson and Bitous [2]; (17) 
and (18) compare with their (25) and (26) and 
show that the characteristics of our partial 
differential equation are reaction paths; (22) 
divided by da, dT gives their equation (82). 
The important difference between our equation 
and theirs is that the sign of dS dr is negative. 
This is because we are integrating from the end 
of the reactor not from the beginning and carries 
with it the advantage that we can determine the 
initial value of a instead of finding it by trial 
and error. To do this we observe that when 
S = 0, f(x, Yo. 0) = 0 and hence f, and f, are 
zero everywhere on this manifold. It follows from 
(18) that the initial value of a is given by 


A (a) = nF (x9)/G (yo) (24) 


But this is just the so-called zero function which 
Amunpson and BiLous use to tell them where to 
stop the integration. Equation (12) immediately 
shows what points are worth considering as 
starting points for 


Ss = naF (x) — b (a) G (yo) 


when S = 0 and a is given by (24); we are only 
interested in regions of the x, yy plane that make 
Ss positive. 

Up to this point we have merely confirmed 
Amunpson and BiLous’ analysis by an entirely 
different method, but this has no great advantage 
unless it be that the theory of first order partial 
differential equations is better known than that 
of functionals. We have still a need for trial and 
error for, whilst they had no means of knowing 
what the initial temperature should be, we have 
no idea of what the values of x, and yy may be at 
the end of the integration. However in the 
matter of restrictions the method of dynamic 
programming has a distinct advantage over 
functional analysis for with the latter it would 
be very hard to say what should be done with such 
an inequality as T, < T < T*. Moreover it is 
essential to impose a restriction of this kind, 
for if G(y = 0, the optimum temperature 
required by (13) becomes infinite ; such is the case 
with a feed of pure A. 

Since a(T) is a monotonic function of T the 
restriction on 7' implies a similar restriction on a. 


(25) 


a* <a (26) 


As soon as the value of a demanded by equation 
(13) goes outside these limits a must be fixed at 
that limit and equation (12) is then integrated 
with constant a (equal to a, or a*). The charac- 
teristic equations for this are again (17)-(21), 
but here a is fixed at a, or a* and the equations 
are easily integrable. In the initial manifold 
S = 0 where a is given as a function of xy yy and 
by (24) we can define two curves I, and I* 
along which the value a is a, or a*. Between 
these curves the initial value of f, is given by 
eliminating a between (24) and (25). Outside 
these curves f, will still be given by (25) but with a 
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fixed at a, or a*. Two other curves A, and A* 
can be drawn along which f, = 0 for a — a, and 
a* respectively. If (b a — db da) does not change 
sign (and this is usually the case) one of the 
curves will lie between the two [° curves. This is 
illustrated in Fig. 1 where A* lies between I’, 
and I* and in this case f, is positive in the whole 
area between A* and the axis y, — 0. Thus the 
region of initial points of interest is between A, 
and the xy axis. 


Fig. 1. The regions of interest on the initial manifold. 


We can now describe the way in which a 
complete solution of our problem is to be generated. 


Let the three regions of the initial manifold be 
denoted by 


I: a=a,, fy>0 
Il: ag<a<a*, f,>0 
III: 


Consider first a point vy = € Y = » in the region 
Il. The initial values for the integration of the 
characteristic equations (17)-(22) are a = &, 
S=f,=f, and a is given by 
A (a) = nF These six equations are 
to be integrated until either a = a* or ay or 
Yo = 0. In case a = a* we continue to integrate 


a=a*, fy > 0. 


the equations (17)-(21) but with a kept constant 
until ry or yg = 0. If (€, 4) be in the region I we 
must integrate the equations with a = a* con- 
stant but substitute the values of 2», yo. f, and f, 
into equation (13) until a = a*, or until x or 
Yq ~ 0 whichever may happen first. In the latter 
case the complete trajectory has been found, 
but in the former case we must continue the 
integration of the full characteristic equations 
(17)-(22) until either a = a* or x or yy = 0, 
and if a reaches a* switch again to the equations 
with constant a = a*. Similarly from region 
III the equations with constant a= a* are 
integrated but with an eye on equation (13) lest 
the unconstrained maximum should fall below a*. 
This however is unlikely for B(a) does not 
normally change sign. Thus a whole region of 
. Yo S space is swept out and for each point 
of this region the value of f(a. yp. 8) and the 
appropriate inlet temperature a (ap, is 


known. If we are given the holding time S and 
the inlet concentrations 2, y, of A and B we find 
the maximum yield immediately as f (2, 8). 
The temperature policy is given by the values 
of a along the characteristic from the point 


Tq: Yo. S to the plane S — 0. For by (17) and (18) 
the characteristic is a reaction path and the value 
of a at a distance s along it is the optimum inlet 
value of a to a reactor of length S — s when the 
concentrations are x (s), y(s) and by the principle 
of optimality this is just the optimum temperature 
at this point. 


3. Consecutive First Orper Reactions 

The foregoing example has been discussed in 
some detail as being sufficiently simple to 
introduce the method yet general enough to be 
interesting. We will now give some definite 
results on the case that Amunpson and BiLous 
discussed before treating more general problems. 
In doing so we can exhibit a device that is often 
useful in reducing the number of variables and 
dispose of two questions which the earlier work 
left unresolved. The results of the last paragraph 
would be two functions f(x, yy, 8) and 
Yo. S) the optimum yield and initial 
temperature of a reactor of holding time S and 
with inlet concentrations x, and yy, and such 


= bo 


functions can be represented as contour surfaces 
in [3] space. This, however, is not convenient 
and it would be of value to work with functions 
of only two variables. In the case of first order 
reactions F = 2, G(y) =y this is possible 
for the equations are homogeneous in the concen- 
trations and we can take z = y 2, their ratio as 
the dependent variable. 

Letting n = 1, F = x, G = y in equations (4) 


and (5) we have 
= (27) 
ar — by, (28) 
whence 
(29) 


dt 


with = = % = Y % when ¢ = 0. The maximum 
yield of B will be proportional to x) so we write 


(Xo Yor 8) = 8 (2, S) (30) 
and since fy = = & — Sy = We 
have in place of (9) the equation 


= max [(a — bz) (1 + g.) — a(g — 2g,)]- (31) 


A common form of the rate coefficients a and b 
isa =k, exp(—E, RT)andb =k, exp(— FE, RT) 
in which case 


b =x (32) 


where r= E, and bb, *” 
In this case A (a) = db/da =(p r)a''”" and 


B(a) = —ra*, and the condition for an un- 
restricted maximum in equation (31) is 
(1 + Be) 
Substituting this into equation (31) gives the 
non-linear equation 
r\*{1 + (1 + 2) 

where R =r (1 —r). The characteristic equa- 
tions are 


a‘ l-r)/r (33) 


(34) 


a{1+(1 
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dS 


(36) 
a(g—1) + pal (37) 
T 
pa" (1 +g), (38) 
a 8s, (39) 
and (40) 
dr 


Formidable though these may seem it is actually 
possible to find a solution in finite terms, for 
dividing (35) by (40) we would have a single 
ordinary differential equation of Riccati’s form 
for z as a function of a. Letting B(a) = 
(6 — a) ra® the standard method of solution 
gives 


{exp B B(a)} 


— —— (61) 
{exp + fo B(a)da 


for the solution for which z = (g when a = ay. 
However this answer would require considerable 
computation for the next step is to substitute in 
(40) and integrate for time 7 as a function of a. 
The modern digital computer is as well adapted 
to the solution of a set of simultaneous ordinary 
differential equations (35)-(40), as it is to the 
integrations required by (41). 

We have to solve (31) subject to the boundary 
condition g (29, 0) = 0, that is to find the surface 
S) that passes through the z, axis. To find 
the initial values that the dependent variables in 
(35)-(40) must take we observe that this condition 
implies g=g,=—S=—0 and for any = 
gs and a are given by (34) and (33) 


(Zo = (1 —r) ay = (1 —r)(r/p (42) 


It is now clear that some restraints must be put on 
a for this would otherwise require it to be infinite 
when z = 0, R> 0. Let the inequality restric- 
tions (26) be imposed. Then we can say at once 
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that only where the formula (33) gives a value 
of a in the range (a,, a*) do the preceeding 
equations (34)-(40) apply. If the value is outside 
then the equation (31) has to be solved with a 
held constant at either its upper or lower bound. 
This would lead to the same set of characteristic 
equations (35)-(39) but with a — a* or a = a*. 
These however are easily soluble and give 


(2, = 


exp(— aS) — exp(— bS) — 1 — exp( 


(43) 


with a and b (a) at one of these values. Naturally 
the region of interest is bounded by the curve on 
which g, — 0 for beyond this point it is not worth 
reacting. Equation (39) shows that g, = 0 is a 
characteristic, and since it passes through the 
point on the z axis at which g, — 0, namely 
a, = a, b,, it has the equation 


byl dg + (dg — exp — (ay — (44) 


We may now distinguish three distinct cases 
according as r = 1, 0 <r <1 and r > 1 which 
we will consider in turn. 


r= 1. 

This case is somewhat special since FE, — EF. 
and the relative rates of the two reactions are 
unchanged by temperature. It follows that if 
the concentrations favour the production of B 
then the greatest possible temperature should be 
used to achieve this as fast as possible. As soon 
as the concentrations lead to the production of 
C at the expense of B the reaction should stop. 
This may be seen analytically from equation (31) 
with b = pa for 


&s = a* {(1 + g, — g) — (p — 1g. + p)} 


whenever it is positive. Then the solution is given 
by (43) with a = a*, b = pa*. 


O<r<l. 


In this case we observe that R — r (1 — r) 
is positive and ranges from 0 to infinity. It follows 
that on the initial line, the 2,-axis, there will be 


three regions : 


I <2, =(r/p)(a*y'/* 
where ay a*, 


< 2* = (r/p)(a,)/" 
where a, < dy) < a*, 


=< % 


<% <Z=(1/p) (ay) 
where dy = a*. 


Fic. 2. Integration of the characteristic equations, 


These are shown in Fig. 2 in the top section ot the 
figure. The corresponding initial value of gy is 
given by either (31) or (42) and is shown in the 
second section. 


I. (gs)o = (1 — (rp 
(gs)o = — 2 
and it can be shown that (g,)y has a continuous 
derivative with respect to 2 . 

The way the integration of the characteristic 
equations proceeds can be illustrated by consider- 


ing three typical characteristic paths. In the first, 
AB, we start from a point A in the region I 


(46) 


= 

—j 
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where the required temperature a) is maximum. 
The initial value of g, is given by (46) and the 
integration of equations (35)-(39) proceeds with 
a = a* until z, becomes zero (the other possibility 
of the optimum a falling below a* does not occur). 
With a characteristic, such as CDE, emanating 
from the region II we start in the region where 
the optimum a is unconstrained and the full set 
of equations (35)-(40) applies with initial values 
given by (45) and (46). The integration proceeds 
until a point D is reached at which a = a* after 
which the characteristic equations for constant 
a = a* must be used until z=, = 0 at EF. When a 
characteristic starts from the third region, the 
first part of it FG is found by integrating the 
equations with a = a,. However a point G will 
be reached at which the criterion (33) for an 
unrestricted optimum a demands a value greater 
than a,. From this point the full equations must 
be integrated until H where a — a* and from 


10 30% 


8 20% 


05 0 20 


Fic. 3. The optimum yield surface. 


then on a is kept constant and equal to “i 
a* until z = 0 at J. The locus of the 
point G at which the change over of § 03 
policy occurs can be found analytically 
for the part of the surface originating 5 62 


from III is given by (43) with a = ay. 
Finally the whole region of interest is 
bounded by the characteristic gy = 0 . 
whose equation is (44). 

The results of an actual calculation are 
shown in Figs. 3 and 4 for the case 

= 0.5, p= 5, a, 01, at 04. In 


14 — 
0 7 
8 

6 

~ ~ 
2 
ie, 


Fic, 4. The optimum temperature policy. (Full lines 
are of constant a; broken lines show some 
charateristics). 


Fig. 3 the contours of the function g (z, 8) are 
shown. From this the maximum yield from 
a tubular reactor with inlet concentrations 


‘Xo Yo and holding time S can immediately 


be obtained as In Fig. 4 the 
characteristics and lines of constant temperature 
policy are shown. From this second figure the 
policy for any conditions x, Yo and S may be 
found. We know that the characteristics are 
reaction paths so that the characteristic from 
Yo to the point where it reaches the 
axis S — 0 gives the relation between z and S$ 
(i.e. between = and ¢) down the reactor. The 
value of a at (2%. 8) gives the optimum inlet 
temperature and by the principle of optimality 
the optimum temperature profile must be given 
by the values of a along the characteristic. 


4 i2 10 8 6 4 2 
Holding time to exit of bed, S 


Fic. 5. Temperature profiles for various lengths of 


25 


reactor and z = 0. 
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Fig. 5 shows some profiles for z, — 0, that is a 
feed free of B. The reactor should not in any 
case have a holding time of more than 13-8 for 
this would put us in the region of negative gy. 
With S = 13-8 the optimum profile is to hold the 
temperature constant at its lowest possible value 
and the yield is 50 per cent. If the reactor is 
shorter, say of holding time 12, it ts necessary 
first to have a very short section of maximum 
temperature followed by a rapidly falling gradient 
and the rest of the reactor at minimum tempera- 
ture; the yield is 48 per cent. With shorter 
reactors the initial section at maximum tempera- 
ture becomes longer and the falling gradient less 
steep ; the yield falls off accordingly. 

Because the absolute optimum for any given 
z is to have a reaction path on the boundary 
Zs = 0 the optimum policy is to have a = a, 
throughout and a reactor of holding time given 
by (44). In particular for a feed of pure A, the 
greatest possible yield is 


be / (ay 
a, 


(47) 


This disposes of a paradox suggested by AMUN Dson, 
that if E, > EF, an infinitely long reactor at zero 
temperature would secure the total conversion 
of A to B. This is indeed the limit as a, +0 
and b, = pa,’ with r < 1, but since tempera- 
tures must necessarily be bounded away from 
absolute zero, equation (47) gives the best that 
can be done with any particular lower bound a,. 


r>1. 

In this case where E, > E, the algebraic 
manipulations of the preceeding case yield the 
same equations. However, we must go back to 
equation (31) to see what they really mean. 
This equation might be written 


fs = Max (aP — 


where P=1+4(1+%)g,—g and Q= px 
(1 + g,). Now in case r > 1 the expression on 
the right hand side, which starts by being zero 
when a = 0, decreases to a minimum and then 
increases monotonically with a. The formal 
manipulations find not a maximum, but this 
minimum so the equations of the last paragraph 


give the minimum yield of B. The maximum is to 
be found by giving a its greatest permitted value 
a* provided this makes gy positive. The known 
solution (43) applies with a — a* throughout the 
region in which gy > 0; this region is given by 
the equation (44) but with a* in place of a,. 
It thus appears that the doubts that Amunpson 
and BiLous express as to the validity of their 
method when r > 1, ({2] p. 126), were in fact well 
founded. 


4 Tue Genera Prosiem or an 
Orrmcum Temperrature Gradient 
We are now in a position to give a general 
solution of the problem of finding an optimum 
temperature gradient in a tubular reactor with a 
number of simultaneous reactions. If n reactions 
are taking place among m (> 2) chemical species 
it is always possible from stoicheiometric con- 
siderations to choose n species and express the 
concentration of the remainder in terms of these 
n and the initial concentrations. Let the concen- 
trations of these species be denoted by ¢,, 
i= 1,2,***n, then the kinetic expressions can 
also be written in terms of these concentrations. 
A balance over a short section of the reactor 
then gives 


de, 


(48) 


with the inlet concentrations specified by 
(0) = (49) 


Instead of maximizing the yield of any 
particular substance we will consider quite a 
general function V (e,, ¢,) which will represent 
the value of the product stream when the con- 
centrations are ¢,,*** ¢,. The problem is to find 
a temperature profile T(t), 0 <¢ <8, such that 
V (e, (8S), eg (8), is as large as possible. 
It is again clear that when such a profile has been 
found the resulting maximum of V is a function 
of the eg and S. Let 


* (S)] = Max V fe, (8). (8)] (50) 
then by the principle of optimality 
SF 8) = Max f(e, (8). (8). 8 — 8) (51) 


— 
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for all s < S, where the maximization is over the 
choice of T(t) in O <t<-s. By letting s +0 
we arrive at the partial differential equation 


ly ) 
= Max R, (ej. 
38 i (Cio 
the maximization being over the possible choices 
of T,, the inlet temperature. In keeping with a 
conventional notation for partial differential 
equations of the first order, let us write 


yf 
deo = Prov 


Then the equation (52) is 


(52) 


Cio = Ss = 


n 
1 
T, being always chosen to maximize the second 
term. If 7, is not restricted it must satisfy 
n 
Ri (a, = 0 (54) 
1 
where R; denotes dR, dT. 
(54) requires a T, outside the acceptable range 
T,. T* then 7, must be held fixed at one or 
other of these values. Moreover it is necessary 
to check that the solution of (54) does correspond 
to a maximum and not to a minimum. If we 
denote )R; da, by the characteristic equa- 
tions of (53) and (54) can be written 


— = — R,(2, T,), i = (55) 
dr 
dS _ 
= p,R, (2,7), (57) 
T 
Ep Ryle. Ti (58) 
jul 
_ (59) 
dr 


and if 7, is not held fixed we must add a further 
equation 


zz p;(R, Ri, — R; Ry, Ry (60) 
j= 


T 


If the solution of 


27 


where as before primes denote partial differentia- 
tion with respect to T. In general these equations 
must be solved simultaneously but we notice 
that if n = 2 the ratio p, p, can be eliminated 
between equations (54) and (60). In this case it 
is only necessary to solve three equations, (55) 
for i = 1 and 2 and the modified (60), 

The initial values for the integration of these 
equations come from the condition that 


(eyo ** * ©) = V * * * Cno) (61) 


It follows that for a given initial point ey = €, 
say, we have for + = 0 


a, = a4, =0 (62) 

é=1,2,°°° 

and p,., and 7, are given by 
» 
Pass = 2 (Se) To) (63) 

0= 


provided that the solution of (64) lies in the 
acceptable range T,, T*. If it lies below 7* then 
T, = T* and p,,, is given by (63) with T, = 7%, 
if above T* then T, = T* and we put 7, = 7T* 
in (63) to give p,.,. Again the value of p,4, = fs 
is constant on a characteristic so the region of 
interest on the initial manifold is bounded by a 
hypersurface of [x — 1] which can be found from 
(63) and (64) by setting p,,, = 0. Within this 
region there are three subregions according as the 
initial value of T, is T,, lies in the range T,, T* 
or is T*. The integration proceedure is the same 
as in the previous cases. If it is necessary to 
start out with T, fixed at one of the bounds only 
the equations (55)-(58) need be integrated, but 
constant watch is kept on the criterion (54) to see 
when the optimum 7, leaves the bound, The 
integration then continues with the full set of 
equations (55)-(60) until 7 reaches the other 
bound or one of the 2; becomes zero, 


5. Discussion 
Whilst the preceeding section gives a method of 
finding the optimum temperature gradient in quite 
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a general case there are other questions that 
come to mind, If for example we include the cost of 
constructing and maintaining a reactor of holding 
time S, then a new profit function should 
be used which takes this into account. Suppose 
that this cost is a function C (S) then we might 
take a profit function P = V[e, (8). ¢, (8)] 
— C(S) and seek to maximize this. Here we are 
going to choose not only T(t), 0 — S, but 
also S itself and the optimum when these choices 
are made will be a function only of eyo. a9. ** * Cpe 
say (Cy This problem can be solved 
by the method of the previous paragraph. All 
we need do is to keep track of 
dP av 

in the integration along the characteristic lines 
and stop as soon as this derivatives reaches zero. 
The value of S at which this step is made will be 
a function of the 1, there, say 
Sa (Cy Cae). Then g (tie = S (Cre 
Cy: 7) — G (a). 

We have not said anything as to the possibility 
of actually attaining a given profile. If @ is the 
rate of heat removal per unit length from the 
tubular reactor and ¢, and p the specific heat 
and density of the reaction mixture a_ heat 
balance gives 


dT 
2 MR eT) 


(65) 
where the H, are suitably defined heats of reac- 
tion. From the required profile Q may be caleu- 
lated as a function of ¢ from this equation. In 
view of the rather steep gradients that may be 
demanded for optimum working the values of 
Q may become too large for normal engineering 
practice to handle. It would therefore be of value 
to find an optimum design for which Q should not 
exceed some maximum, Q < Q*. This can be 
done by varying the upper limit of permitted 
temperature 7*. The maximum value ot Q will 
be a function of T*, say Q,, (7'*), and this may be 
determined quite readily for several 7T*. The 
solution is then given by choosing the T* that 
makes Q,, (7'*) = Q*. Alternatively the value of 
Q at each point of a given trajectory may be 


calculated from (65) during the course of the 
integration. Then as soon as Q reaches the value 
Q* the temperature T must not exceed the value 
it has then reached. By suddenly requiring the 
temperature to be constant, i.e. the left hand side 
of (65) to be zero, a discontinuity in Q is produced 
which lowers it below the maximum permitted 
value Q*. 

If a tubular reactor with a definite cooling 
system, such as co- or counter-current tube 
cooling, is to be designed for optimum yield by 
choosing the correct inlet temperature and area 
of cooling tube we are actually faced with a 
much more difficult problem than choosing the 
whole temperature profile. This is because the 
problem is too determinate, and fixing the inlet 
temperature and the cooling tube area fixes the 
whole temperature profile. It is hoped to discuss 
this type of problem elsewhere. 

Finally it should be evident that other features 
of tubular reactors can be considered by the 
methods of dynamic programming. With gas 
reactions the pressure will be an important 
variable and we would work with partial pressures 
rather than concentrations. With catalytic 
reactions there is also the possibility of varying 
the grade of catalyst and its selectivity throughout 
the bed. This problem can also be treated by 
dynamic programming. 

Acknowledgements—1 am indebted to Professor N. R. 
Amunpson for many valuable discussions of this subject, 


and to Mr. D. Rupp of this University for calculations 
on which Figs. 3 and 4 are based. 


NOTATION 
Ala) = db/da 
a = rate constant of A > B =k, exp(— E,/RT) 
Bia) = (b — Aa)/ A’ 
b = rate constant of B-> C =k, exp(— E,/RT) 
¢; = concentration i = 
C(S) = cost of reactor of length S 
F(x) = dependence of rate A —> B on concentration 
of A 


S(%q Yq S) = optimum yield from reactor of holding time 


G(y) = dependence of rate of B — C on concentra- 
tion of B 
k, = constants in reaction rates a and b 
n = number of reactions 
P = profit 
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Q = rate of heat removal Vic) = value of product 
r= E,/E, a = concentration of A 
Re=r/(l —r) y = concentration of B 
R, = rate of reaction z=a/y 
S = total holding time of reactor Bla) = (b — a)/ra® 
t = current holding time along tubular reactor p = kak," 
T = temperature 7 = variable along characteristics 


(1) 
[2] 
(3) 
(4) 
[5] 
[6] 
[7] 
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Abstract—A solution is presented for the diffusion of heat or mass within a developing laminar 
boundary layer in which the velocity is described in terms of directional components, cach raised 
to an arbitrary power. Reduction of the general solution to specific cases is demonstrated. The 
Schmidt or Prandt! number exponent must be between values of 4 and 4 depending on the ratio 
of the concentration to hydrodynamic boundary layer thickness. 


Résumé—Les auteurs étudient la diffusion de chaleur ou de mass a linteriéur d'une couche limite VOL. 


laminaire en formation; la vitesse est exprimée en fonction des composantes spatiales affectées 
dune puissance arbitraire, [ls montrent comment passer de la solution générale aux cas particu- 
liers. Les exposants des nombres de Prandtl et de Schmidt doivent étre compris entre § et 4 
suivant le rapport de lépaisseur de la couche limite de concentration ou de température, a 
lépaisseur de la couche limite hydrodynamique. 


Zusammenfassung— Es wird cine Ldsung angegeben fiir Stoff- oder Wirmediffusion innerhalb 
einer sich entwickelnden laminaren Grenzschicht, in der die Geschwindigkeit durch Richtungs- 
komponenten dargestellt wird, die durch je eine tibergeordnete Kraft hervorgerufen werden. 
Die Anwendung der allgemeinen Losung auf spezielle Fille wird aufgezeigt. Der Exponent der 
Schmidt- oder Prandt-Zahl, der von dem Verhiltnis der Konzentration zur Dicke der hydro- 
dynamischen Grenzschicht abhingt, muss sich dabei zwischen | und } bewegen. 


Whereas the process of diffusion within a The following simplifications are now imposed ; 
laminar boundary layer has been analysed 
extensively for cases involving various velocity te ow 
distributions [1, 2], the general instance in which (2) Y negligible contribution of the terms 
the velocity at the fluid-solid boundary is (de/ dy) and 


described in terms of 2 and y raised to arbitrary phe error involved in neglecting the perpendicular 
powers, } and « respectively, would appear to be component of velocity will be discussed later. 
worthy of analysis. Thus the velocity relative 

to that at the boundary-layer edge is assumed : 

to be of the form c(O,y) =e, 


a) b (a, 0) = Cy 
where the variables « and y are being expressed 
as Reynolds numbers in the direction of flow The equation to be solved is then 
and perpendicular thereto. 
The diffusive process is considered to be one (° = 
of heat or mass transfer which may be described > v¥} oF 
by the equation X 


de de Ye Ye 
dy * [Fs (2) 


(1) Constancy of physical properties. 


Y = 


SOLUTION 
Let z = X* Y* 
Therefore A(z) = ¢(X, ¥Y) = A(X" Y"). 


Performing the necessary differentiation and 
substituting into equation (3) yields 
dA 
K xe 
de 


dA va we 
(= E (n — 1) dz (4) 
m= —b—1 


equation (4) reduces to 


aA 
—— = (a 


(5) 


where 


B= —(n—1)/n 


and P = B = fluid properties group 


Letting p = dA dz and integrating, one obtains 
dA dz = 2 

The constant C is obtained by further integration 

over the entire range of the variable z. Thus 


Recognizing that the integral on the right-hand 
side of equation (6) is of the form of the gamma 
function, I’, one finds that 


AA (— 2)**! 
Therefore 
Tat) 


The net diffusive flux at the boundary, y = 0, 
then becomes, in terms of # and y, 
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def 
U_ 
(- ) (8) 
Since — == kAe 
y=0 


then the instantaneous value of the transfer 
coellicient, k, is the bracketed expression in 
equation (8). 

The average transfer coeflicient may be ex- 
pressed in terms of the exponents a and 6, the 
Schmidt or Prandtl number, P = v D, and the 
length Reynolds number, Re = U, x v, as 


[K (b + 1)/(a + (a + 
U (Py a~—1)/(a+2) (Re)! b—1)/(a+ 2) (9) 


Discussion 
A number of cases pertinent to the science of 
chemical engineering following by specification 
of the exponents, a and 6. 


Case 1 
The velocity is invariant in w and y, that is 
a=b=Oand K=1 
+ OD 2/2 
(1/2) \4(a@ U) 70 
which is the Hicere [3] relation for the average 


transfer coetlicient, where «/U = @, the time of 
transient diffusion. 


k (10) 


Case 2 
The velocity is a simple linear function of y 
and independent of x, that is 


b=0, K =(v U,) 


@ = i, 


This specification will be recognized as the 
Leveque [4] approximation for heat transfer 
between the wall of a conduit of diameter 2yo, 
and a fluid in laminar flow, where the average 
velocity is U,. In this instance equation (9) 


reduces to 
k — 1-615 (5) (11) 
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In terms of the Nusselt number, 


2.1/3 
Nu — k D — 1-615 (“see (12) 


Equation (12) is identical to that derived by 
Leveaue [4] and by and Proroxp [5}. 


Case 3 

The velocity is a function of both « and y, 
in a fashion dictated by laminar boundary layer 
theory. That is 


a=i, 6=1/3 


and the velocity distribution becomes 
ur 


UU, Ky Ji 


The term AK, will be recognized as the first, 
linear term of the various power series often 
employed to describe the velocity distribution 
within a laminar boundary layer on a flat plate. 
Equation (9) becomes, in this case, 


k=0617U_ P** 


Ky (13) 


(14) 


Equation (14) has been successfully employed 
in the prediction of fluid-particle mass transfer 
in fixed beds (6), by defining the distance over 
which the boundary layer develops in terms of 
the particle diameter. 

The error involved in ignoring the term 
V(Qde dy) in the above development may be 
assessed by a comparison of equation (14) with 
the result found when this term is included. 
Following Scnuicurine [2] 


4 


(15) 


=1—2e + 
Ca 

where 

%q 

5, 8 = and = 5-83 (16) 


€ = 


Differentiating (15) to obtain the flux at y = 0, 
and solving for k in a manner comparable to that 


outlined earlier, one obtains 
k = 0-687 UL Ret (17) 


Hence, the error involved by the assumption 
underlying equation (14) is about 20 per cent. 


Scumipr or Pranpt. NumbBer 


FUNCTIONALITY 

Equation (9) states that the exponent on the 
Schmidt or Prandtl] number must lie between 
4 and 4. It seems reasonable to suppose that 
Case 1 above, in which the velocity is invariant 
in « and y, corresponds to that circumstance in 
which the major resistance to the transport 
process lies outside the hydrodynamic boundary 
layer. This behaviour may be expected when 
the diffusivity is rather high (low Schmidt or 
Prandtl number). On the other hand Case 3 
corresponds to the instance wherein the major 
resistance lies within that portion of the boundary 
layer which can be represented by the linear 
approximation stated in equation (13), 

It is of interest to inquire into the exact 
nature of the fluid property number dependency 
as a function of the ratio of the concentration or 
thermal to the hydrodynamic boundary layer 
thicknesses. 

This ratio is detined as 


A = 3,/8, 


where according to ScuLicwrine [2], A is specified 
in terms of the property group P as 


== 


313 P (18) 


and 

H(A) = 554 ino > + Teo 4 for A <1 
and 


H(A) 


10 


21 
15 

1 1 
180 

(The third term in the latter equation is 
erroneously printed as 9/15(1/A*) in the 
SCHLICHTING text.) 

A plot of the solution to the above equations 
is shown in Fig. 1, which indicates that the pro- 
perty number exponent tends toward a limiting 
value of } at values of that group of less than 
1-0. Jn terms of k, 


k= O817U, Re? P> 10 
k= 126U, P< 0-05 
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Diffusion within a developing boundary layer — a mathematical solution for arbitrary velocity distrib:ution 


Fic. 1. Ratio of concentration to hydrodynamic 
boundary layer thickness as a function of Schmidt 
number, 


It may be concluded that in the property value 
regime of from 0-05 to 1-0 no unique exponent 
upon the Schmidt or Prandtl number prevails. 


NOTATION 


Transformed concentration variable 
Exponents in equation (1) 
Concentration of diffusing substance 
Constant defined by equation (7) 
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Perforated-tube flowmeter 
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Abstract—A cheap and simply-constructed visual flow-meter is described. It consists of a 
Rotameter-type bob fitting fairly closely inside a perforated tube. Most of the fluid flows out of 
the holes below the bob; the height of the bob indicates the flow-rate. A simple theoretical 
treatment is given. Following experiments with air and water, formulae are given which can be 
used to specify the leading dimensions of a meter for a given duty. 


Résumé — Les auteurs décrivent un débitmeétre a lecture directe facile réaliser et peu codteux. 
Il consiste en un flotteur du type Rotamétre bien ajusté dans un tube perforé, La plus grande 
partie du fluide s‘écoule par les trous au-dessous du flotteur. La hauteur du flotteur indique 
la vitesse d’écoulement. Suit une étude théorique simple. Aprés des expériences avee lair et 
Neau, les auteurs donnent des formules pouvant étre utilisées pour trouver les dimensions essen- 
tielles dun appareil de mesure pour une plage déterminée. 


Zusammenfassung — Ein billiges Durchflussmessyeritt cinfacher Konstruktion wird beschrieben. 
Es besteht aus cinem Schwimmer vom Rotameter-Typ mit geringem Spiel innerhalb cines per- 
forierten Rohres. Der grésste Teil des strOmenden Mediums fliesst durch die Ldécher unterhalb 
des Schwimmers aus. Die Hohe des Schwimmers zeigt die Durchflussgeschwindigkeit an. Eine 
cinfache Erlauterung der theoretischen Grundlagen wird angegeben. Aufgrund von Experimenten 
mit Luft und Wasser werden Formeln entwickelt, die zur Berechnung der Hauptabmessungen 
eines Messgeriites fiir cine bestimmte Durchflussleistung dienen konnen. 


Tue instrument described here is a simple concentric. The inner tube must be vertical, 
visual flow-meter, which can be constructed from — but this setting is not critical unless the bob is a 
cheap materials without the need for any preci- very close fit. 


sion machining. L 


Description or FLOwMETER 

The flowmeter is shown diagrammatically in 
Fig. 1. It consists of two concentric transparent 
tubes. The inner one is made of plastic (e.g. 
Perspex) and perforated by a regular array of 
holes. A “ bob,” similar to that used in a Rota- 
meter, is suspended by the stream of fluid in the 
central tube. A proportion of the fluid (depending 
on the closeness of fit of the bob) leaks past the 
bob ; the remainder flows out through the holes 
in the central tube. The height of the top of the 
bob indicates the flow-rate of the fluid. 

A cross-wire is fitted to the top of the inner 
tube, to prevent the bob from being blown out, 
and a coiled-spring arrester at the bottom. A 
wire spacer can be fitted to keep the two tubes 


13 
96 

H 

4 

imi 

Fic. 1. Flow meter. 
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Thick-walled Perspex tubes tend to crack after 
drilling. This can be prevented by annealing at 
70 C for a few hours. 

A sphere can be used instead of a Rotameter- 
type bob, but is usually less satisfactory because of 
instability. It would be more satisfactory in 
principle to have a continuous opening in the 
walls of the central tube rather than a series of 
holes. Vertical and spiral slits were tried, but 
neither were satisfactory. It seems that non- 
precision methods of fabrication do not produce 
slits of uniform width and finish, and that the slit 
may widen, temporarily or progressively, if the 
inner tube bulges under pressure. 

The minimum flow-rate which can be measured 
is that which leaks past the bob ; this can be made 
smaller by making the bob a closer fit in the tube. 
The range of the instrument depends on the 
weight of the bob and the number and diameter 
of the holes. If the bob is hollowed-out its weight 
can be varied by putting in fine lead shot. The 
accuracy depends on the size of the holes ; 
other things being equal, twice as many holes, 
each of half the cross-sectional area, will roughly 
halve the average error. 


EXPERIMENTAL 


Experiments were carried out with air and 
water, using a variety of tubes and bobs; in 


most cases the latter were hollow and weighted 
with varying amounts of fine lead shot. Typical 
results are shown in Fig. 2. Table 1 shows the 
range of variables covered. 


50 


20 
4, cm 


Fic, 2. Typical calibration curves. 


Water Air 
N = 63cm"! N =63em"! 
d = 1275 em d = 1275 em 
A = 0-011 em? A = 0-011 em? 
D = 1-19 ¢m = 119 ¢em 
= 2-35 wt = 139g wt 


Table 1. Showing range of variables covered 


T 


W /981 
(Weight of Q 
bob immersed | (Maz, values 
in fluid) measured) 
(g wt.) see) 
017-068 
0-26-0-76 


12 34 
480-770 


0-83-1°56 
1-29-1-82 


0-90-2°55 
1-39-2-18 


1-39-2°55 


Chem. Engng. Sci. Vol. 13, No. 1. August 1960. 


2000 
5 
25 7 woo § 
CY: 
} 
N d | D A Q% 
(em=) (em) (em) (em?) | Fluid | | (cm? sec) 
63 685 055 oo | water | 245 
| | | | air | 120-130 
1-19 0-011 te 27-44 48 
an 
| | 0-00126 | air 
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Comparison Turory 


The rate of leakage, Q,. past the bob will be 
the same (for a given fluid) at all flow-rates. 
Following the conventional treatment for the 
Rotameter, we have 


D 


where D and d are the bob and tube diameters 
respectively, p the density of the fluid and W 
the submerged weight of the bob (foree units). 
K is related to the discharge coetlicient for tlow 
past the bob; for all the experiments discussed 
here the value of A was found to be between 
O-8 and 1-2. 

The non-linearity of the relation between the 
flow-rate, Q, and the height of the bob, MH (see 
Fig. 2) is due mainly to the fact that the fluid in 
the central tube declerates as it flows upward, 
and that its pressure therefore rises. A simplitied 
analysis is given in the Appendix, and leads to 
the relation 


nd? 


where Y — @ — Q, is the total flow-rate through 
the holes beneath the bob, No the number of 
holes per unit height and 4 the area of one hole. 
C is the discharge-coeflicient for the holes below 
the bob, and & a constant relating to pressure- 
recovery due to decleration. This relationship, 
with C = land k — 1, is compared with certain 
of the experimental curves in Fig. 3; the remain- 
ing experimental curves, for cach hole diameter, 
lie between the curves shown. It will be seen that 
the results conform broadly to the * theoretical ” 
expression, equation (2) above. The expression 


NAH 


0-5, k 0-56) fits all the 


(2) 


(corresponding to C 


experimental curves to within a factor of 1:5 
over most of the range. 

Equation (2) implies that at a certain flow-rate, 
given by 


-+ 
Fic. 3. Experimental results. 
N d | D A W981 
Curve Fluid | (em!) (em) | (em) (em?) (g wt.) 
1 water 0-635 0-546 0-675 
2 air 63 0-635 0-537 oon 0-259 
3 water 0-635 0-546 ool 0-170 
4 air 6-3 1-275 1-19 ool | 1-90 
5 air 6-3 1-275 1-14 0-011 1-82 
6 air 6-3 1-275 1-19 0-00126 1-4 
air 6-3 1-275 000126 1-39 
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2D kp 
the bob will be blown right up the tube, no matter 
how long the latter may be. This condition is 
approached by certain of the curves in Fig. 3. 


CONCLUSIONS 
The following expressions may be used to 
caleulate approximately the leading dimensions 
of a perforated-tube flowmeter for a given duty, 
at least so long as conditions do not depart too 
widely from those described, Final adjustments 
are easily made by adding lead shot to the bob. 


(4) 


ep 
Qo = D (5) 


It follows from equation (6) that the relationship 
between Q and H will be almost linear at lower 
tlow-rates (when 2NAH wd? < 0-5, i.e. when the 
area of the holes below the bob is less than the 
cross-section of the tube). 
D (7) 


< wd*/4, 


where 


NAH 


There is an upper limit, Q,,,,. to the flow-rate 
which can be measured, regardless of the length 


of the tube : 
all 
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APPENDIX 


The general problem of flow in uniformly-tapped pipes 
is a difficult one ; a much simplified treatment is used 
here. 

The following approximations are made : 


1. Wall-friction is neglected ; 
2. The discharge coeflicient through each hole below 
the bob is taken to be the same ; 


The velocity below the bob is regarded as negligible 
compared with that at the entrance to the inner 
tube ; 

The bob is regarded as simply a flat disk, its position 
in the tube being that of the top of the real bob ; 
It is assumed that pressure-recovery due to decelera- 
tion in the inner tube can be expressed by 


dP = — kpd (u*) (9) 


where * is a constant of order of magnitude unity [1]. 


The pressure in the space between the inner and 
outer tubes and above the bob is taken to be 
uniformly equal to Po. 


Let the pressure in the inner tube at a height h above 
the level at which the perforations begin be P. Then 
the rate of discharge through each hole at this level 
(from assumption 2) is 


where A is the area of each hole and C the discharge 
coeflicient. Then, if q is the residual vertical flow in the 
inner tube at height A, and N the number of holes per 
unit height, 


dh 


The pressure P’ at height H (just below the bob) must be 
sufficient to support the weight of the bob, so 


(10) 


=W (11) 


where D is the diameter of the bob. 
u in the tube is related to q by 


The linear velocity 


(12) 


Substituting for «u in equation (9), and integrating be- 
tween A (an arbitrary height) and 1/7, we have 


whence P — P,= (13) 


(<4 


Substituting this expression for (P — Py) in equation 
(10), and integrating between O and H, we obtain an 
expression for the total rate of flow, Q’, through the holes : 
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[(8 p) — 2K 


sin AN (14) 


(In setting the upper limit of the right-hand integral 
equal to 0, we are introducing assumption 3). 

This treatment is easily extended to take into account 
the acceleration of the fluid in the annular space between 
the inner and the outer tubes, which is ignored in assump- 
tion 6. The left-hand side of equation (14) should be 
multiplied by a factor which is approximately 


where a, is the cross-sectional area of the inner tube and 
ay that of the annulus. In the experiments described here 
the correction never amounted to more than about 
1 per cent. 


REFERENCE 


NOTATION 


Area of hole 

Cross-sectional area of perforated tube 
Cross-sectional area of annulus between tubes 
Discharge coetlicient for hole (equation 10) 
Interior diameter of tube 

Maximum diameter of bob 

Height above lowest hole 

Height of top of bob above lowest hole 
Constant in equation (1) 

Constant in equation (9) 

Number of holes per unit height 
Pressure 

Pressure outside tube and above bob 
Pressure at height // 

Vertical flow-rate at height 4 

Total tlow-rate 

Flow-rate past bob 

Flow-rate through holes 

Velocity through holes 

Submerged weight of bob (force units) 
Density of tuid 
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Note on a method for computing the effect of temperature on reaction yields 


(Received 24 April 1960) 


Rastoe: and Densicn [1] have pointed out that the effect 
of temperature on the yield of a chemical reaction which 
is not allowed to proceed to equilibrium is quite different 
from the effect of temperature on the equilibrium yield. 
To compute this effect in a particular case they have 
integrated the rate equation which describes the detailed 
kinetics of their example, Subsequent elementary manipu- 
lations yield an expression for r = (dA, dT) (dA, dT) which 
may be as high as 10 for A A, — 0-9 Here A is the degree 
of advancement of the reaction at a fixed time, and A, 
is the degree of advancement at equilibrium. 

The object of this note is to develope an approximation 
which enables one to effect a quite considerable simplitica- 
tion in the calculation of r by circumventing the solution 
of the compileated differential equations which describe 
the reaction kinetics. The subsequent calculations also 
become appreciably simpler and the cost in accuracy is 
slight. 

We write the chemical reaction in question as 


X,, stands for 1 mole of component « and vy, is the stiochio- 

metric coefficient of a. Stoichiometric coeflicients are taken 

as positive for products and negative for reactants, 

Denoting the molar concentration of component « by n,. 
the degree of advancement, A, is defined by 


dn, =v, dA (2) 


The constant of integration implied by equation (2) is 
determined by setting A = 0 at the initiation of the 
reaction, The rate equation is written 


(1) 


dA 
= T) (3) 
where T is the absolute temperature. The form of g 
depends on the particular reaction under investigation. 
The approximation considered here is the utilization of 
the observation that not too far from equilibrium g (n,.T) 
can be expanded about its equilibrium value of zero, and 
only the term linear in A retained [2]. Thus we write 

approximately 
dX 


dt 


~ 


— 


where the subscript ¢ indicates that the derivative is to 
be evaluated for the equilibrium values of the n,. The 
relevant solution of equation (4) is 


Writing f, simple differentiation yields 


din A, /dt 


(¢) 


d\/dT 

dd, /aT 
(dind,/dT) may be determined by thermodynamics, 
(din | /dT) depends on the mechanism of the reaction 
under investigation. 


As an example we consider the reaction also chosen by 
Rasrocr and Densicu 


X, +X, 


(6) 


-X, =0 (7) 


with n, (0) = ng (0) = a, ng (0) 
the rate law is 


= = 0. Assuming 


dX 


= ky (nyng — ngny) (8) 


where K is the equilibrium constant, the relaxation time 
is given by 


= 2ak,K'? (9) 


1 


ak! ? 
= ith) 
re 1+ Ki? om 
Then from equation (6), the van’t Hoff equation and the 
Arrhenius equation, 
) (11) 


where E is the activation energy for the forward reaction, 
and A// is the enthalpy change of the reaction. Numerical 
results are given in Table 1. The results of Rastoer and 
Dennicu obtained by exact integration of the rate equation 
are also tabulated for comparison. Note that we do not 
consider values of f less than 0-9, for our considerations 
yield a reasonable approximation only near equilibrium. 

Table 1 also gives calculations for two other rate laws, 
together with values determined by exact integration. 


Table 1. 
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[r (enact)). 


Comparison of r computed by present method [1 (present)| with r computed by the exact method of Reference 1, 


Reaction 


Kinetic 
type 

E (Weal) 

H (Keal) 


kK 


Bimolecular 


bimolecular 


| 


cis trans 
CH CH | CHCH 
3 3 


unimolecular 

18 
[8] 


CH,CO,C 
CH,COOH 
Unimolecular 

himolecular [4] 


[6] 


0-005 [5] 


208 


r (present r (exact) r (present) r (exact) r (present) r (exact) 


286 2:8 2:8 2-88 


240 


580 re tes is 


Numbers in brackets are literature references. 


treatment of Rasroc: and Denrnicn is about the same as 
in the temperature case. 


In the case of opposing uni-molecular reactions the approx- 
imation we use is in fact exact since the actual rate equa- 
tion is linear to start with. In all cases shown the approx- 
imate calculations vield numbers in tolerable agreement 
with the exact calculation. 

The method presented here can be extended without 
difficulty to cases where several competitive reactions 
occur simultaneously. The differential equations are 
again trivial, but the subsequent algebra is quite compli- 
cated and the results do not appear to be in a particularly 
useful form. In addition, the method can be used to treat 
the effect of other parameters such as pressure on reaction 
vields. Such effects were also discussed by Rasrog: and 
Denpsicn. No new features are introduced by considering Mu 
pressure rather than temperature, so we shall not go into f 
details. The numerical aggreement with the more exact r 


Novation 
- molar concentration of component « 
> stoichiometric coeflicient of component « in a 
reaction 
forward rate constant 
equilibrium constant 
progress variable of a chemical reaction 
relaxation time of a chemical reaction 
= activation energy for the forward reaction 
enthalpy change 
fractional attainment of equilibrium 
ratio of (dA dT) at fixed time to (dA, dT) 
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Book Reviews 


Polymer. Butterworths Scientific Publications, London. 
Published quarterly. Annual subscription £5, $15.00, 


During the last twenty years, with the growth in im- 
portance of the chemistry and technology of polymers, 
the need for an outlet for the large volume of published 
matter on these topics has become clear. The need was 
recognised some fifteen years ago in America with the 
appearance of the Journal of Polymer Science and more 
recently with the addition of the Journal of Applied 
Polymer Science, by which it is hoped to separate the more 
fundamental and applied aspects. “ Polymer,” a new 
journal produced by Butterworths Scientific Publications, 
is designed to cover both these aspects and its choice of 
editorial board should go far towards achieving a sound 
balance. The journal is to appear quarterly and the first 
part of Vol. T comprising some 120 pages, has just appeared. 
Containing twelve full papers on the chemistry and physics 
of a variety of polymers, in solution and otherwise, and 
on polymerization processes, as well as short letters and 
three book reviews, the objectives of the journal are 
satisfactorily reached in this issue, Nor should it be difficult 
to maintain this level from the large amount of relevant 
work which is proceeding in industry, research institutions 
and university laboratories. The concentration of polymer 
work into a single journal will undoubtedly be welcome 
to the polymer worker, but the price of the publication 
(£5 per annum for less than 500 pages) will ensure that it 
does not occur widely in personal libraries. 

P. JoHnson 


Introduction to 
MeGraw- Hill, Lon- 


W. H. Corcoran and W. N. Lacry: 
Chemical Engineering Problems. 
don 1960, 185 pp. 53s. 6d. 


M. G. Lantan: Fundamentals of Chemical Engineer- 
ing Operations. Constable, London 1960, 644° pp. 
62s. Gad. 


TeN YEARS ago it was possible to own most of the text- 
books in Chemical Engineering which existed, and to be 
fairly well acquainted with the others. Nowadays they 
are being published so rapidly that unless complimentary 
or review copies arrive on one’s desk one may never see 
the newcomers, 

Most of these textbooks naturally arise from courses 
taught at various universities, and each therefore reflects 
the character of the syllabus and the intellectual atmos- 
phere at a given university. For this reason, the majority 
are of limited value as aids in teaching elsewhere, Certain 
of the textbooks published before the flood began achieved 
widespread recognition ; for instance, Principles of Chemical 
Engineering by Wacker, Lewis, and GILLILAND 


of M.L.T. ; Elements of Chemical Engineering by Bavorn 
and McCane of Michigan ; Chemical Process Principles 
by Hovucen and Watson of Wisconsin ; and Chemical 
Engineering by Covutson and Ricnarpson of Imperial 
College ; however, later arrivals have proved less generally 
useful. Indeed, present tendencies in the teaching of 
chemical engineering, which stress development from the 
scientific fundamentals rather than description and analysis 
of the “unit operations,” call for monographs on the 
principles and applications of chemical and engineering 
thermodynamics, fluid mechanics, transport processes 
ete., rather than textbooks which purport to cover the 
major part of the subject. 

It is, in fact, misleading to give a textbook a title which 
suggests that it deals exhaustively with chemical engineer- 
ing (unless it confines itself to a superficial survey, or 
consists of many volumes). Chemical engineering, to my 
mind, includes chemical thermodynamics and_ kinetics, 
fluid mechanics and transfer processes, plant design and a 
number of other subjects. If a book leaves out some of 
these subjects its title should reflect the fact. 

The book by Professors Corcoran and Lacey consists 
of chapters on units, engineering calculations and measure- 
ments, behaviour of gases (without reference to fluid 
mechanics), material and energy balances, introductions 
to chemical equilibria and kinetics, and the manufacture 
of ammonia, nitric acid, sulphuric acid and sodium 
hydroxide. It is concerned principally with the applica- 
tions of elementary 
engineering problems, and not at all with such subjects 
as fluid mechanics and transfer processes. It would not 
be possible to guess this from the title. 

Probably because of the particular arrangement of the 
svilabus at the California Institute of Technology, the book 
is designed “for use before physical chemistry is 
encountered.’ There is probably more to be said for this 
approach to engineering education than would generally 
be admitted in English universities. The principles of 
thermodynamies, for instance, are much more difficult 
than the applications; should we perhaps start with 
the simple formulae of thermodynamics and their applica- 
tions, before returning to the first principles of the subject 
and developing it in all its rigour? This method would at 
any rate re-enact the historical development of thermody- 
namics, and at the same time make it more digestible to 
many. Reaction kinetics, as the authors rightly remark, 
is still (at any rate, from the chemical engineer's viewpoint ) 
an empirical subject, and loses very little by being dealt 
with as such. 

In English circumstances this book might well prove 
useful when a first-year introductory course in chemical 
engineering is being taught in parallel with a course in 
physical chemistry. It is well written, and contains 


physical chemistry to chemical 
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over 50 useful examples. Unfortunately, it contains one 
or two fallacious or misleading statements, particularly 
with regard to the kinetic basis of equilibrium and the 
use of activities instead of concentrations in rate equations, 
which can only store up trouble unless the student is 
warned against them. 

Professor Larian’s book contains no chemistry ; he 
might claim that the word “ operations in the title 
should warn us to expect none, Even so, in attempting 
to cover fluid mechanics, transfer processes and the 
“unit operations ” in some 600 pages, he has shown that 
one volume is inadequate for the purpose. This would 
have been a more valuable book, in my opinion, if all 
fundamentals had been taken as read, and the space 
available devoted to their application to the “ unit 
operations,” and to the description of the principles of 
operation of various types of equipment. As it is, the 
treatment of the fundamental principles (e.g. of diffusion) 
is perfunctory, The author emphasises in his preface that 
his book has a practical bias. The syllabus at Michigan 
State University presumably reflects his contention that 
the graduating engineer “ should possess the minimum 
practical knowledge to qualify for his first job.’ This is 
an aim with which one cannot quarrel, although it must 
be admitted that the graduate who is trained to be 
immediately useful will probably not also have time at the 
university to master the principles of the subject with 
sufficient certainty to play a leading part in radically new 
technological developments. However, there is certainly 
room in industry for graduates of both kinds; the im- 
portant thing is that the “ practical” course should 
be as sound, in its own way, as the ~ intellectual ~ 
course. Judged by this criterion, the weakness of Professor 
Larian’s book lies in a lack of precision and clarity in his 
discussion of the fundamentals of such subjects as fluid 
mechanics and transfer processes. There is no need to be 
imprecise merely because the discussion remains at the 
elementary level ; in fact it is just at this level that it is 
important to be crystal clear. 

Professor Larian is much more successful in his 
discussion of the general features and the behaviour of 


various types of equipment, and in his presentation of the 
conventional methods of calculation ; from this point 
of view the book may prove to be a useful text for the 
more “ practical type of chemical engineering course, 
although one would like to see more attention paid to the 
question of how far the conventional methods can be 
expected to give the right answer in practice. 


P. V. 


F. Rumrorp: Chemical Engineering Materials. 
(Second Edition) Constable, London 1960, Octavo, 384 pp., 
Gd. 


Tur cnorcr of a material of construction for a chemical 
plant remains for the most part a matter for experience, 
and this experience differs in suprising ways from one 
industry to another. Theory about corrosion helps little, 
and the body of practical knowledge is vast and not easy 
to classify. It is, therefore, a pleasure to welcome the 
second edition of the excellent work on this subject by 
Dr. Rumrorp, who manages to arrange a wealth of 
information in a small volume and in such a way that 
it is easy to discover the properties of a specific material. 

The text is divided according to the materials used for 
construction, and the principal mechanical properties 
and the resistance to corrosion of each material is discussed. 
Information is also given on constructional methods and 
and costs. Metals and inorganic materials, common and 
uncommon, are treated, and there are preliminary chapters 
on the mechanism of corrosion and on the testing of 
materials, (The book by Ler, Materials of Construction 
for the Chemical Process Industries, which is classified 
according to the chemicals it is desired to contain within 
the constructional material is complementary to Rumford’s 
book. The two together form a good library on this 
subject). 

In this second edition the author has included more 
information on the cost of material than in the first 
edition, and there is a different treatment of the high 
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